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ABSTRACT 

The dust ejecta of Main-Belt Comet P/2010 A2 (LINEAR) have been ob- 
served with several telescopes at the at the Observatorio del Roque de los Mucha- 
chos on La Palma, Spain. Application of an inverse dust tail Monte Carlo method 
to the images of the dust ejecta from the object indicates that a sustained, likely 
water-ice driven, activity over some eight months is the mechanism responsible 
for the formation of the observed tail. The total amount of dust released is esti- 
mated to be 5x10^ kg, which represents about 0.3% of the nucleus mass. While 
the event could have been triggered by a collision, this cannot be decided from 
the currently available data. 

Subject headings: minor planets, asteroids: general — comets: general — comets: 
individual(P/2010 A2 (LINEAR)) — methods: data analysis 

1. Introduction 

On January 6, 2010, P/2010 A2 (LINEAR), a comet-like object, was discovered by the 
LINEAR sky survey. The orbital elements of the object (a = 2.29 AU, e = 0.12 and i = 5.26 
deg) are typical of an inner Main-Belt asteroid belonging to the Flora collisional family, and 
suggest that it is unlikely to have originated in the classical comet source regions (i.e., the 
Kuiper Belt or the Oort cloud). The object reached perihelion, at a heliocentric distance of 
r=2.0 AU, on 4 December 2009. 



Observations taken a few days after its discovery showed an inactive nucle us lying outside 
a dust cloud that looks like a cometary tail without central condensation (ILicandro et al. 
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2010al lbl: | Jewitt et al.ll2010l ). This suggested that the observed dust cloud and the nucleus 
were the result of an impact between two previously unknown asteroids: the dust cloud is a 
plume of dust and the nucleus is what remains from the largest of the asteroids that collided. 

Owing to its cometary-like aspect and orbital pa rameters, P/2010 A2 can be classified 
as a Main-Belt Comet (MBC, iHsieh fc JewittI l2006l ). In contrast with the other known 
MBCs, however, P/2010 A2 has a significantly smaller semi-major axis (2.29 AU versus 2.7 
AU and ~3.2 AU). Until now, the activity observed in MBCs was found to be compatible 
with a water-ice driven activation mechanism suggesting that those asteroids retained ice 
layers below their surface and, under certain conditions, become "activated asteroids" . The 
presenc e of water-ice on the surface of 24 Themis, t he parent asteroid of the Themis family 
MBCs fjCampins et al.ll2010l : iRivkin &: Emeryll2010l ) strongly support this. 



One of the activation mechanisms could be a collision between two asteroids. Collisions 
are known to take place regularly, but they are so rare that none of the dust plumes that 
they should generate has ever been seen. If P/2010 A2 is the debris of a collisional event, this 
would be the first time that the ejecta from such a collision is observed soon after it happened. 
It would provide a unique opportunity to learn something about asteroid collision processes 
and about the internal composition of asteroids. Alternatively, if the observed activity is 
sustained in time as in comets, this leaves open two interesting problems: (1) if water-ice 
sublimation is the activation mechanism, how does water-ice survive in an asteroid with such 
a small semi-major axis; or (2) is there any other mechanism capable of ejecting dust in a 
similar manner as water-ice sublimation does ? 

In this paper we present and analyze images of P/2010 A2 obtained with three telescopes 
at the Roque de los Muchachos Observatory (ORM), La Palma, Spain. In Sect. 2 the 
observations and data reduction are presented. In Sect. 3 we analyze the combined images 
using the inverse Monte Carlo dust tail fitting method (e.g. iMoreno et al.l 120041 : iMoreno 
20091 ) to study possible ejection scenarios. The conclusions are presented in Sect. 4. 



2. Observations and data reduction 



Cepa et al. 


2000; 


Cepa 


2010) 



camera-spectrograph at the Gran Telescopio Canarias (GTC), with the Auxiliary camera- 
spectrograph (AC AM) at the William Herschel Telescope (WHT), and with the Andalucia 
Faint Object Spectrograph and Camera (ALFOSC) at the Nordic Optical Telescope (NOT) 
all located at the ORM. 
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The observational circumstances are given in Table 1. For the observations with the 
NOT, the telescope was set at the comet's rate motion, while sidereal tracking was used for 
both WHT and GTC In all cases the images were obtained in service mode by telescope 
staff and were reduced by subtracting the overscan level and flat fielded using standard 
procedures. 

The observations with OSIRIS at the GTC were made with the Sloan g',r', and i' filters. 
OSIRIS provides a field of view 7.8' x 7.8' with a gap of 9.2" in the middle and a pixel scale of 
0.125"pixel~^. To increase the signal-to- noise the data were binned in 2x2 pixels. The images 
were calibrated using photometric zeropoints determined from standard star observations. 

The observations with ACAM at the WHT were made with the Sloan g' and r' filters. 
ACAM is mounted permanently at a folded- Cassegrain focus of the telescope and has a 
circular field of view of 8' diameter with a pixel scale of 0.25 "pixel"^. The resulting calibrated 
WHT images show brightness levels in agreement with those obtained at GTC within 10%. 

The observations with ALFOSC at the NOT were made with standard Johnson- Cousins 
V and R filters. ALFOSC provides a field of view of 6.5'x6.5' with a pixel scale of 0.19 
"pixel"^. These images were calibrated with stars in the field of view using magnitudes from 



the USNO-Bl.O catalog (IMonet et al.ll2003l ). which provides a photometric accuracy of ~0.3 
mag. 



Using the transformation equations of iFugukita et al. (Il996[) and the magnitude of the 



Sun in the standard Johnson- Cousins V filter (l0=-26.75. ICoxll2000l ). we derive rQ=-26.96. 
Since i?0=-27.29, the r' images of the object obtained at GTC and WHT were transformed 
to R standard Johnson- Cousins magnitudes by subtracting 0.33 mag, where we assume for 
the object the same spectral dependence as for the Sun within the bandpasses of these two 
red filters. The resulting GTC, WHT, and NOT R magnitudes are consistent with their 
errors and we conclude that the object did not change brightness significantly during our 
observations. 

The available images in each night were combined in order to improve the signal-to- 
noise ratio and converted to solar disk intensity units appropriate for the analysis in terms 
of dust tail models. Figure 1 depicts some representative images obtained on the different 
dates with the instruments mentioned above. In this figure one can clearly see the inactive 
nucleus (marked with an arrow) lying outside the dust cloud that looks like a cometary tail. 
The absence of a dust cloud surrounding the nucleus indicates that the dust emission has 
stopped before the observations. Using aperture photometry, we determined the nucleus R 
magnitudes for the GTC, WHT, and NOT images whenever possible (see Table 1). The 
aperture size was of 1.5", and the nucleus profiles were stellar-like. For the GTC images. 



- 5 - 



which have the best S/N ratio, and using the formahsm by lBowell et al.l (119891 ) with a slope 
parameter of G'=0.15, we obtained a nucleus absolute magnitude of iJ=21.3±0.3. Hence, a 
nucleus diameter of D=220±40 m can be determined, assuming a bulk albedo of p=0.11, 
typical of a S-type asteroid. This assumption is based on the fact t hat the S-type aster oids 
are the most common objects in the inner Main Asteroid Belt (e.g. iBus fc Binzelll2002l ). 

The combined images from GTC and WHT obtained on the 17th and 21st January 
show the best S/N ratio (see Fig. 1), and were selected for a more detailed analysis. To 
carry out the analysis in terms of dust tail models an additional rebinning was made of the 
images. The WHT image was rebinned to 4 pixels, giving a spatial resolution of 772.4 km 
pixel"^, while the GTC image was rebinned to 3 pixels to give a resolution of 576.3 km 
pixel"^. Finally, the rebinned images were rotated to the {N,M) coordinate system, where 
M is the extended radius vector from the Sun, and N is perpendicular t o M and directed 
opposite to the object's motion along its orbit (IFinson fc Probsteinlll968l ). 



3. Application of the inverse dust tail model to P/2010 A2 



We have performed an analy sis of both the WHT and GTC ima ges by the inverse Monte 

Moreno et al.l l20oi: iMorend 120091) . which is based on the 



Carlo dust tail fitting code (e.g. 



Fullelll989l . |2004| ). Briefly, a large amount 



same procedures as that developed by Fulle (e.g. 
of particles (typically 10^) of a selected size range are released from the surface of the object 
in a given time interval, with an assumed velocity law, which may be a function of time 
and particle size. The particles are then submitted to the radiation pressure and gravity 
field of the Sun, and their Keplerian orbits are computed. The ratio of the force exerted by 
the solar radiation pressure and the solar gravity is given by /3, which can be expressed as 
/3 = CprQpr/{2pr), where Cpr.=1.19x 10~^ kg m~^, and p is the par ticle density. For p article 
radii r > 0.25 pm, the radiation pressure coefficient is Qpr ~ 1 (IBurns et al.lll979l ). The 
inverse method involves the inversion of the overdetermined system of equations AF = I, 
where A is the kernel matrix containing the dust tail model, i.e., the surface density of the 
sample of particles integrated over time and (3, F is the output vector which contains the 
time-dependent distribution of (3, and I is the brightness of the input image in the selected 
region of the photographic (A^, M)-plane. 

The physical appearance of the tail, detached from the nucleus and sun-ward oriented 
(see Figures 1 and 2), suggests that most of the dust was already ejected sometime well 
before the date of the observation, and should be composed of relatively large particles, 
otherwise they would have been blown away by radiation pressure. This is clearly demon- 
strated in Figure 2, where the combined WHT image is shown along with the synchro nes 
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( iFinson &: Probsteirull968l ) corresponding to ejection times from 272 days before perihelion 
to 3 days after perihelion. Assuming a particle density of 1000 kg m~^, these synchrones 
would correspond to particle sizes between 0.001 cm and 1 cm. To start the simulations, 
and based on the synchrone map, we assumed initially that dust has been ejected at any 
time starting approximately a year before the observations, with sizes in the range 0.001 cm 
to 1 cm. We fixed the particle density to p= 1000 kg m~^, and the albedo times the phase 
function to be Ap=0.04, i.e., a Halley-like value. In principle, we assumed isotropic ejection. 
The better spatial resolution of the GTC image compared to the WHT image might suggest 
anisotropic ejection, which is particularly indicated by the shape of the isophotes near the 
nucleus. We will come back to this point later. 

For lack of better information, we assumed an ejection velocity as given by f oc 
which from hydrodynamical considerations is typical from gas drag by sublimation processes. 
We did not assume any temporal variation in the ejection velocities, which is reasonable 
taking into account the low eccentricity of the orbit. We found that the best fits to the 
images are found when the velocity is given hy v = 1100/3^/^ cm s~^. This im ply a velocity 



of 8.5 cm s~^ for a grain of r=l cm. Using the formulae by IWhippld (119511 ) for a nucleus 
of 220±40 m in diameter as given above and heating efficiency factor of 0.1 we get, at a 
heliocentric distance of 2 AU, a range of velocities of 12-14 cm s~^ (and smaller values for 
lower values of the heating efficiency), which are not far from our model results. This would 
then be consistent with a typical cometary activity scenario where the outgassing is driven 
by water-ice sublimation. 

The result of the model fit to the WHT image is shown in figure 3(a), and the derived 
dust mass loss rates as a function of time are shown in figure 4. The model isophotes 
reproduce quite closely those of the WHT image. The variation of the mass loss rate indicates 
that the onset of the activity occurred in late March 2009, with a maximum of activity around 
early June 2009, with a peak of dust loss rate of about 5 kg s~^, and decreasing to 0.1 kg s~^ 
near perihelion in early December 2009. The integrated dust loss is 5.4x10^ kg. Assuming 
a spherically-shaped nucleus of diameter of 220 m, and a bulk density of 3000 kg m~^, the 
dust mass in the tail corresponds to 0.3% of the total nucleus mass. 

We tried other model combinations to search for other possible solutions. In particular, 
we tried to fit the data by assuming narrower time intervals for the dust ejection, which could 
support an interpretation of a single collision event without any further outgassing activity, 
but we could not find any acceptable solutions even modifying also the particle sizes and 
ejection velocities with values well above an d below the expected values that occur in the 



laboratory experiments such as presented by iNakamura fc Fujiwaral ( ll99ll ). 



If the isotropic ejection model with the same input parameters is applied to the GTC 
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image, we obtain the result shown in figure 3(b). In this case, due to the fact that a 
certain portion of the tail is clipped, the dust mass loss differs from that found for the 
WHT image, particularly in that the mass loss becomes zero at times earlier than 180 days 
pre-perihelion (see figure 4). For the same reason, the integrated ejected dust mass is also 
lower at 1.7x10^ kg. Looking closely at the fitted isophotes, some mismatch in the tail 
head near the nucleus can be seen, and also in the outermost isophote towards the sun-ward 
region. Therefore, we decided to incorporate ejection anisotropy in the model by considering 
emission from a selected active area at a particular region of a spherically symmetric rotating 
nucleus. The rotation period was set to 3 hours. After trying tens of models with varying 
rotational parameters and location of the active area, a very good fit was obtained by taken 
the argument of the sub-solar meridian at perihelion as $=180°, and the obliquity as 1=95°, 
with an active area occupying a part of the southern hemisphere of the object in the latitude 
range [-90°,-40°]. The resulting fit to the GTC observation is shown in figure 3(c). If a 
collision triggered the outgassing event, it could have occurred at high southern latitudes on 
the object. Although this model seems highly speculative, and other alternative scenarios 
exist, it does produce an excellent fit to the data. 

For the resulting time-averaged size distribution, a fit indicates a power-law with an 
index of -3 .4±0.3, a value typically found in the analysis of dust ejected from comets (e.g. 



PuUdlSQM). 



The above fits indicate that the observations cannot be explained by dust released in 
a sudden event such as a collision and subsequent ejection of asteroid surface material, but 
it does need a sustained cometary-like activity during several months after the onset of the 
activity as described above. Whether the onset of the activity was triggered by a collision 
event or by another mechanism that might be similar to those occurring in other Main-Belt 



Comets (IHsieh fc Jewittll2006l ). can not be decided from the available data. 



4. Conclusions 

The application of the inverse dust tail Monte Carlo fitting method to images of the tail 
of P/2010 A2 tail has revealed that the observed brightness distribution is only consistent 
with a sustained (possibly water-ice driven) cometary-like activity that spanned approxi- 
mately 8 months, which is estimated to have started in late March 2009. The total amount 
of the ejected mass was about 5x10^ kg, when assuming an albedo times phase function of 
4%, and a particle density of 1000 kg m~^. It is not possible to asses whether the outgassing 
event was triggered by a collision or by another unspecified mechanism. However, the de- 
tection of a sustained cometary-like activity over a considerable period of time has obvious 
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consequences for the internal structure of the object. It might indicate that some subsurface 
ice layer exists in this object. Another option is that there exists another mechanism which 
is capable of causing ejection of dust in a manner similar to water-ice sublimation. 

A more sophisticated (but also more speculative) dust tail model indicates that the 
dust ejection was localized and originated on the southern hemisphere of the nucleus, with 
rotation parameters $=180°, obliquity /=95°, and a rotation period of ~3 hours. These 
parameters can be confirmed with further follow-up observations. 
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P/2010 A2 (UNEAR) 




JAN 14 (NOT) 




JAN 17 (GTC) 



JAN 21 (WHT) 



JAN 23 (NOT) 



Fig. 1. — Images of P/2010 A2 (LINEAR) obtained at different dates with different telescopes 
at the Observatorio del Roque de los Muchachos in La Palma, as described in the text. The 
telescope abbreviations are NOT = Nordic Optical Telescope, WHT = William Herschcl 
Telescope, and GTC = Gran Telescopio Canarias. The images have been rotated in order to 
have the tail along the X-axis. The Sun direction is approximately toward -|-X. Whenever 
the nucleus is visible, it is marked with an arrow. 
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N [km] N [km] 

Fig. 2.— Left: Resulting calibrated and rebinned P/2010 A2 WHT image on the {N,M) 
plane, acquired on January 21st, 2010. On the right, the synchrone map is overplotted. 
Synchrones are shown at (in counterclockwise sense) +3, -52, -107, -162, -217, and -272 
days since perihelion date, with maximum ^=0.0595 (corresponding to a minimum particle 
radius of 0.001 cm for a density of 1000 kg m""^). 
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Fig. 3. — Results of the dust model applied to the WHT and GTC images. Coordinate axes 
correspond to the (A^, M) system (see text). Black contours, with isophote levels of 5xl0~^^, 
10"^^, and 2xl0~^^ solar disk intensity units, correspond to the observations. Red contours 
correspond to the model. Panel (a): The isotropic model applied to the WHT image. The 
physical dimensions of the images are 57154 km x 194632 km. Panel (b): The isotropic 
model applied to the GTC image. The physical dimensions of the images are 25932 km x 
48407 km. Panel (c): The anisotropic model applied to the GTC image, with same physical 
dimensions as those of panel (b). 
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Fig. 4. — The derived dust mass loss rates versus time from the WHT and GTC images 
using isotropic and anisotropic dust ejection models as indicated. The arrow marks the date 
of the observation of the WHT image. 
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Table 1. Observational circumstances of P/2010 A2. Position data include heliocentric 
distance r, geocentric distance A, and phase angle a. Last column shows the 
Johnson-Cousins /^-magnitude of the nucleus with their estimated uncertainties. 



UT date Telescope r (AU) A (AU) a° Nucleus i?-mag 



14 Jan 2010 


NOT 2.56 


m 


2.013 


1.039 


5.53 


23.3±0.4 


17 Jan 2010 


GTC 10.4 


m 


2.014 


1.047 


7.20 


23.5±0.3 


21 Jan 2010 


WHT 4.2 


m 


2.016 


1.061 


9.32 


23.2±0.3 


23 Jan 2010 


NOT 2.56 


m 


2.016 


1.065 


9.87 





